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Abstract 

We study the kinetics of phase transitions in a Rayleigh-Benard system after 
onset of convection using 2D Swift-Hohenberg equation. An initially uniform state 
evolves to one whose ground state is spatially periodic. We confirmed previous 
results which showed that dynamical scaling occurs at medium quench (e = 0.25) 
with scaling exponents 1/5 and 1/4 under zero noise and finite noise respectively. 
We find logarithmic scaling behavior for a deep quench (e = 0.75) at zero noise. 
A simple method is devised to measure the proxy of domain wall length. We find 
that the energy and domain wall length exhibit scaling behavior with the same 
exponent. For e = 0.25, the scaling exponents are 1/4 and 0.3 at zero and finite 
noise respectively. 
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1 Introduction 


The dynamics of pattern formation in systems away from equilibrium is a fascinating 
subject to study. There are numerous examples that exist in nature, such as Rayleigh- 
Benard convection, chemical reaction and biological pattern formation etc. Much work 
|jl|] has been done on Rayleigh-Benard convection due to its relative simplicity. There, a 
fluid is conhned between two horizontal plates which are heated from below. When the 
temperature difference is big enough such that the Rayleigh number R exceeds a critical 
value Rc, the fluid in the uniform state becomes unstable to one which consists of spatially 
periodic convection rolls. These rolls form domains with a typical size which grows with 
time. Swift and Hohenberg developed a simple model p to describe this process. The 
order parameter equation reads. 
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eb — ~ (1 + + V 


where b is related to the vertically averaged magnitude of the velocity component normal 
to the plate. The reduced Rayleigh number e = {R — Rc)/Rc measures how far the 
system is above the onset of convection. The thermal noise r] satishes the usual relation 
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{ri{x,t)7]{x',t')) = F6{x — x')5{t — t'). This model captures the three basic features of 
pattern formation in the Rayleigh-Benard system, namely, initial growth of the instability, 
its nonlinear saturation and development of a spatially periodic pattern. 

Away from the threshold where e ~ 0(1), one can perturb around the ground state 
which is periodic in x and define a phase 0 as in 

V' ='0o[0(x,t)] +0(0, 

where V0(x, f) = q(x, f) and the gradients of local wavenumber q are of the order of a 
small number 0 Thus one obtains the phase equation [^,0. To lowest order in e, 

dt(j) = D±d^(j) + O||c}jj0. 

where T and || denotes directions along and normal to the rolls and the O’s are diffusion 
coefficients. Direct dimensional analysis tells us that the characteristic length grows as 
However, D± —> 0 due to local wavenumber relaxation. The next order expansion 
gives us a fourth order gradient term and one expects a scaling behavior L ^1/4 

in 

the transverse direction (along the rolls) Bi which physically corresponds to curvature 
relaxation |^. 

Scaling behavior is seen in numerical simulations. However, the exponent is 1/5 for 
noiseless case and is apparently changed to 1/4 when finite noise is added Non¬ 
potential variants of the Swift-Hohenberg system which take into account of mean held 
how give the same result B (non-uniform horizontal motion of the huid, important for 
huid with small viscosity). 

The fact that the phase equation predicts a diherent value of exponent from simulation 
and that the exponent is noise dependent is intriguing. Actually, since the phase equation 
starts from perturbation around an ideal state, it fails to encompass the many defects 
that are present in a real system. However, the coarsening process in pattern formation is 
brought about by defect movement. Through defect annihilation, larger domains of nearly 
ideal conhguration form at the expense of smaller ones. This important role played by 
defects is well-known; a good example is the kinetics of the formation of the Meissner 
state in a type H superconductor at zero held B, where the inter-vortex potential U{r) 
determines the scaling behavior of the inter-vortex distance. 

2 Numerical Results 

We approach the problem primarily as numerical experimentalists. We try to measure 
strategic quantities which could shed light on the mechanism of the pattern coarsening, 
especially quantities which could be related to defects in the system. We use an explicit, 
hrst order accurate, pseudospectral method. The boundary is periodic to simulate large 
aspect ratio systems. Two e values are used. For e = 0.25, the lattice is 512 by 512, time 
step size if 0.1. For e = 0.75, the lattice is 256 by 256, time step size is 0.05. In either 
case. There are 8 lattice points per ideal period. In the case of e = 0.25, we also did a case 
with hnite noise to study the effect of noise on the dynamics. To obtain the statistics, 
each run is repeated with 10 different initial random configurations to average over. 


A. e = 0.25. 

We measure three quantities. One is the structure factor S{k). This is the canonical 
measure from which one can calculate the dynamical scaling exponents. We also measure 
the energy of the system as well as the sum of all the domain wall lengths. The former 
could signal, in this relaxational model, whether bulk or defects are contributing most to 
energy relaxation, while the latter could give us information on defect dynamics. 

The structure factor is defined as the statistical average of the fourier spectrum of the 
equal time correlation function, averaged over all directions of the 2D wave vector. 


S{k) = FT ■ (V'(x,t)V'(x',f'))- 

The Swift-Hohenberg system belongs to the instability class |^. Accordingly, S{k) 
is peaked at k = ko away from the origin, corresponding to the ground state and most 
unstable mode above onset of convection. Due to the nonlinear interaction, there are peaks 
corresponding to higher harmonics present in S{k). Away from these higher harmonics, 
it is observed that S{k) is a function of — /cg|. We £t S{k) to squared Lorentzian form 
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S{k) 


S{ko) 

(i + (r(P-fc2)2)2) 


where S'(fco) is the peak amplitude and T is a fitting parameter from which one can 
calculate L, the characteristic width of S{k). The two parameters S'(fco) and L exhibit 
scaling behavior. The exponents are consistent with results obtained by earlier workers. 


L ~ S{ko) ~ where 0L = | ^ ^ (1) 

where we define the discretized noise = F/{Ax)‘^, where Ax = vr/d is mesh size used 
in our simulation. The data plots for S'(/co) and L are shown in Fig 1. 




Fig 1. L/S(ko) vs Time ( = 0 left and = 0.05 right). Left and right reference 

lines have slopes 1/5 and 1/4 respectively. 





Using these characteristics, we scale the structure factor from different times. The 
central region of S{k) collapses onto a universal curve (Fig 2). 
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Fig 2. Scaled S(k) vs - kQ. Noiseless. 

It is interesting to notice that S'(fco) scales as L [|] not as is the usual case for phase 
ordering in a 2D system. This is due to the fact that, in the process of Fourier trans¬ 
formation S{r)^S{k), most of the contribution comes from annulus k^dk instead of area 
kdk at the origin where k ~ dk. Since F is basically L, we write the scaling form of S{k) 
as 


S{k) = LI ■ f{Ll\e - kl\) 

where I is a constant length to give the right dimension, proportional to the periodicity. 

There is an energy functional H in our model, where 

//= i|(l + 4)^1= - 

The system evolves in a way such that the energy is decreasing as a function of time. The 
questions then are: is there energy scaling? If so, what is the scaling exponent? And is 
it related to the scaling of the characteristic length L? Scaling in energy was observed in 
other systems such as phase ordering in nematic films [^] where it is found to be mostly 
due to bulk energy relaxation. 

The ground state energy density can be easily calculated to be — to the lowest 
order in e. In the noiseless case, this is also the long time limit. We can measure the 
excess energy E of the system relative to that of its ground state. For the noisy case, 
it is a bit tricky, because there is a ‘thermal’ energy component from the external noise 
source. So the long time limit is unknown. What we did was to determine this value 
empirically by fitting the energy evolution to a power law allowing for an offset. Then 
excess energy is found by subtracting this offset. Whether this is a good fit is answered 
by seeing whether the scaling prevails over several decades. 

Indeed, a nice scaling regime is found. However, the scaling exponent is not related 
to that of the characteristic length L in any simple way (such as L~^ as expected from 



ordinary bulk relaxation). We think the dominant contribution comes from defects as will 
be discussed below. Our results are shown in Fig 3, and may be summarized as: 


E t where 


1/4 , Ed = 0, 

0.3 , Ed = 0.05. 
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3 Energy (left axis) 


3 DWL (right axis) 
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Fig 3. E/DWL vs Time (F^ = 0 left and = 0.05 right). Left and right reference 
lines have slopes 1/4 and 0.3 respectively. 


An accurate measurement of domain wall lengths should be very helpful in clarifying 
one of our main concerns: what is the role of defects? However, it is quite difficult to 
determine where a domain wall is located mathematically. Other researchers have used 

H. 


optical hltering procedures to visualize domain walls 


However, no attempt has been 


made to measure the length of these domain walls. Furthermore, the filtering procedure 
is very computationally intensive. 

We devised a simple method to identify domain walls and obtained a proxy of the 
domain wall lengths. The method is based on the observation that, within a domain, the 
order parameter configuration is close to the ground state solution, 'ip ~ Asin(ko • x). 
We have VpJ ~ Ako cos(ko • x) and ~ + (V'^)^//co- The converted amplitude A is 

almost constant everywhere (Fig 4) except at grain boundaries, which can thus be 







located by filtering with respect to A. We set up a threshold so that, if the calculated 
A'^ is bigger than 0.7 x Max{A^) + 0.3 x AVG{A^) or smaller than 0.7 x Min{A^) + 0.3 x 
AVG{A‘^), that point is counted as belonging to a domain wall. The length proxy DWL 
is dehned as a count of all such points. 

Scaling behavior is observed for DWL. The exponent is different from that of the 
characteristic length L. Instead, it is the same as that of the energy relaxation (Fig 
?)){DWL for noiseless case was an average of 5 runs since we devised the above way of 
measuring DWL after the hrst 5 runs). This strongly suggests that excess energy is 
primarily distributed at domain walls and other defects. This part of the energy relaxes 
slowly, via defect annihilation. On the other hand, bulk energy is the fast mode and 
scales as ^^d decays much faster than the excess energy E. At this point, it 

is tempting to conclude that defects are indeed the driving force behind the coarsening 
process due to its dominant contribution to the excess energy. However, it is still not 
clear how one could link the evolution of E and DWL with the growth of L. 

DWL ~ where c^dwl = | Z g 05 (3) 


B. e = 0.75. 

For large e, the system is much further away from equilibrium. We anticipated that 
there would be more isolated dislocations and they could give rise to a different coarsening 
process. This is indeed what we found. The system evolves slowly. The characteristic 
length exhibits logarithmic scaling behavior, while the excess energy decays even slower 
and the number of defects stays almost constant (Fig 5). 

L ~ Siko) ~ logt, Ed = 0. 
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Fig 5. L/S(ko) vs Time (left) and E/DWL vs Time (right) 

A heuristic argument for the logarithmic behavior uses the fact that the pair poten- 
between two dislocations is exponentially decaying with respect to the distance 


between them, namely. 







Setting r X z, one has, for the damped defect motion, 


dr dH 

^ dt dr 




( 5 ) 


The leading order of the solution to above equation gives us logarithmic behavior. Now 
substitute this leading order behavior r ~ log(t) into the defect pair potential energy, one 
hnds that the pair potential energy part scales as much faster than what one sees in 
the E vs Time plot shown above. Thus, there must be other signihcant energy sources 
which could be some energy background not important to dynamics. 


3 Conclusion 

In this paper, we presented results from our numerical simulation of the Swift-Hohenberg 
system. For a medium quench (e = 0.25), the system has many domain wall defects. 
We conhrmed the dynamic scaling behavior in structure factor observed by previous 
researchers. We showed that the energy of the system also exhibits scaling behavior. 
Its scaling exponent does not have a simple relationship with that of the characteristic 
width, making it difficult to draw any concrete conclusion. However, one notices that the 
energy relaxation is slower than what one would expect from bulk contribution, namely, 
L~^. We devised a simple way to measure the proxy of domain wall length. The DWL 
has the same scaling exponents as the energy, suggesting that energy is concentrated on 
defects. This is consistent with the notion that defects are the driving force behind the 
coarsening process. We hud that adding noise into the system speeds up the evolution, 
because all the exponent values increase. 

For a deep quench (e = 0.75), we see a much slower coarsening process where there 
are mostly isolated dislocations in the system. The characteristic length scale shows 
logarithmic behavior. To explain this, we give a heuristic argument assuming overdamped 
isotropic dynamics under the influence of the pair potential between defects. 
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